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Abstract. An enhanced sulfate aerosol layer has been observed near 25 km 
accompanying springtime ozone depletion in the Antarctic stratosphere. We use a one- 
dimensional aerosol model that includes photochemistry, particle nucleation, 
condensational growth, coagulation, and sedimentation to study the origin of the layer. 
Annual cycles of sunlight, temperature, and ozone are incorporated into the model. 
Our results indicate that binary homogeneous nucleation leads to the formation of very 
small droplets of sulfuric acid and water under conditions of low temperature and 
production of H 2 SO 4 following polar sunrise. Photodissociation of carbonyl sulfide 
(OCS) alone, however, cannot provide sufficient S0 2 to create the observed 
condensation nuclei (CN) layer. When subsidence of S0 2 from very high altitudes in 
the polar night vortex is incorporated into the model, the CN layer is reasonably 
reproduced. The model predictions, based on the subsidence in polar vortex, agree 
with in situ measurements of particle concentration, vertical distribution, and 
persistence during polar spring. 


Introduction 

Balloon-borne counters have been able to measure aero- 
sols with radii greater than 0.01 ^m in the stratosphere since 
1973 [Hofmann and Rosen , 1977; Rosen and Hofmann , 
1977]. These fine particles are generally referred to as 
condensation nuclei (CN) because of their small size and 
their role as sites for condensation of sulfuric acid vapor 
[Hofmann and Rosen , 1985; Hofmann , 1990]. Measure- 
ments over Antarctica have revealed an enhanced CN layer 
between 20 and 30 km with a maximum particle concentra- 
tion exceeding 100 cm' 3 at about 25 km [Hofmann and 
Rosen , 1985; Hofmann , 1988; Wilson et al., 1989]. The layer 
has been observed to appear repeatedly in each spring from 
September to November. Hofmann et al . [1989] found that 
the CN layer and the severe ozone depletion (ozone hole) in 
the south polar region are concurrent, and both CN concen- 
trations and ozone abundance are sensitive to local temper- 
ature. Moreover, coincident with the CN layer and the 
ozone hole is an observed decrease in concentration of 
“optical-active aerosols” (OA) with particle radii greater 
then 0.15 pm. Recent measurements in the Arctic during 
winter show a similar layering effect [Hofmann , 1990]. 
Stratospheric CN, which are most likely liquid sulfuric acid 
and water solutions, are believed to serve as sites for 
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heterogeneous chemical reactions. Thus it is of great interest 
to investigate the mechanisms which lead to the formation of 
these CN layers. 

Several mechanisms for the formation of the CN layer 
have been proposed. Rosen and Hofmann [1983] suggested 
that sulfur-bearing gases (carbonyl sulfide (OCS)) could 
accumulate in the polar stratosphere during winter due to the 
lack of sunlight and produce high concentrations of H 2 S0 4 
vapor in spring after sunrise. However, Oppenheimer [1987] 
argued that photodecomposition of OCS is a slow process 
because of the large zenith angle in the polar region; the 
production of H2SO4 through OCS photolysis is not suffi- 
cient for CN formation. Thus he proposed that photooxida- 
tion of OCS or reaction of OCS with elevated levels of OH 
might account for an enhancement in sulfuric acid vapor 
leading to the formation of the CN layer. 

Because CN layers are always observed to be accompa- 
nied by low temperatures, Hofmann [1988] suggested that 
low temperatures in the winter polar stratosphere support 
homogeneous nucleation and rapid growth of fine particles 
due to excessive supersaturation of H 2 S0 4 vapor. Hamill et 
al. [1990] studied the homogeneous nucleation of ambient 
sulfuric acid at polar winter low temperatures and concluded 
that even with modest acid concentrations the formation of 
H 2 S0 4 /H 2 0 droplets would be efficient. However, such acid 
concentrations may not exist in polar winter because pro- 
duction of H 2 S0 4 through photochemical processes ceases 
after polar sunset, while the acid vapor still condenses to 
preexisting particles as temperatures are decreasing. Hamill 
et al. did not include the loss of H 2 S0 4 to preexisting 
particles. Their model did not reproduce the observed sea- 
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Table 1. Photochemical Reactions of Sulfur Components in the Stratosphere 



Reaction 

Rate 


/?> 

Kx 

S + Oj — ► SO + 0 

K { = 2.3 x 10" 12 



Kz 

SO + O; — * S0 2 ■+■ O 

K. = 2.6 x 10“ 13 e 

- 2400/7 

*3 

*3 

SO : + OH + M — HS0 3 + M 

K o ,0 ° = 3.0 x 10 -31 , 
K™ = 1.5 x 10‘ 12 , 

, n = 3.3* 
, m = 0* 

* 4 

K* 

HS0 3 + 0 2 — S0 3 4 - HO, 

AT 4 = 1.3 x 10 -12 e 

- 330/7 

*5 

K5 

S0 3 + h 2 o -+ h 2 so 4 

K 5 = 6.0 x 10’ l5 t 



*6 

OCS 4 O — » SO 4 - CO 

K b = 2.1 x 10“" e 

- 2200/7 

J i 

Jl 

OCS 4 hv — * S + CO (A < 228 nm) 



h 

Ji 

S0 2 +hv -* SO 4 O (A < 220 nm) 




*Three-bo<ly reactions with K 3 = (AT 0 [A/]/(I + tfoEW**)) x 0.6 (1 where 

K 0 = K^°{m00)- n and K x = Kl 00 ( 77300) . 

tThis reaction may require a surface or a small ion cluster to balance energy and momentum. 
Nevertheless, the reaction is apparently very rapid under atmospheric conditions. 


sonal nor altitudinal characteristics of the CN layers. Rosen 
and Hofmann [1983] proposed that evaporation of sulfate 
aerosol in a stratospheric warming event followed by suffi- 
ciently rapid cooling could be responsible for the CN forma- 
tion. However, the evaporation of aerosols may not release 
sufficient acid vapor for CN formation because the concen- 
trations of preexisting aerosol near 30 km are very low. 

Measurements of long-lived tracers suggest organized 
downward motion within the polar night vortex [Toon et at 
1989; Russell et at 1993]. For example, data from the 
Halogen Occultation Experiment (HALOE) instrument on 
the Upper Atmosphere Research Satellite (UARS) show that 
methane (CH 4 ) mixing ratios in polar stratosphere are as low 
as those in the mesosphere. Moreover, the data suggest that 
air with lower CH 4 mixing ratios in the polar vortex does not 
mix with the stratospheric air surrounding the vortex until 
the vortex breaks down in spring [Russell et at ., 1993]. It can 
be inferred that the air at the top of the springtime polar 
stratosphere originates from the mesosphere by descending 
motion in the polar vortex associated with radiative cooling 
during the winter. S0 2 may achieve a high mixing ratio in the 
mesosphere through photolysis of H 2 S0 4 . Although accu- 
rate measurements of the H 2 S0 4 ultraviolet absorption 
spectrum are not available, photodissociation at wavelength 
below -250 nm is likely. S0 2 generated in the mesosphere 
may then be transported into the stratospheric polar vortex 
by subsidence. There, S0 2 has a relatively long lifetime 
because of the absence of OH during polar winter. In this 
situation the mesosphere can be considered as a sulfur 
source for CN formation in spring after sunrise. Hofmann et 
at. [1988] also suggested that a possible source of sulfuric 
acid may be transport from sunlit regions in lower latitude 
with entrance into the vortex above 30 km. Then subsidence 
and radiative cooling in the lower stratosphere may result in 
particle formation. 

Dust particles and micrometeorites from the ablation of 
micrometeors in the upper atmosphere may contribute to the 
CN concentrations above 20 km [Hunten et al. , 1980; Turco 
et al 1981]. These particles can affect the properties of 
stratospheric aerosols by serving as condensation nuclei for 


sulfuric acid droplets, enhancing the sulfate aerosol concen- 
trations and reducing average particle sizes when sulfuric 
acid becomes supersaturated. 

Polar stratospheric clouds (PSCs) may also modulate CN 
populations. PSCs can scavenge preexisting sulfur particles 
and thus decrease the total aerosol surface area in spring. 
The decreased aerosol surface area limits H 2 S0 4 vapor loss 
onto the preexisting particle surfaces and maintains high 
H 2 S0 4 vapor concentrations for homogeneous nucleation 
and CN formation. 

In this study we elaborate on these various hypotheses by 
conducting a series of numerical experiments that include in 
a consistent way chemistry, aerosol microphysics, annual 
variations of temperature and sunlight, micrometeorite 
fluxes, and polar vortex dynamics. 

Model Description 

To simulate the behavior of CN in the springtime polar 
stratosphere, we incorporated sulfur photochemistry and mi- 
crophysics in a one-dimensional version of the tracer model 
developed by Toon et al. [1988]. Table 1 gives the sulfur 
photochemical reactions which are important in the lower 
stratosphere. Under volcanically quiescent stratospheric con- 
ditions, the sulfur oxides originate primarily from OCS trans- 
ported through the tropopause, as first suggested by Crutzen 
[1976]. After reaching the stratosphere, photodecomposition of 
OCS (J j) releases sulfur atoms (S), which are quickly con- 
verted into sulfur dioxide (S0 2 ) through R } and /? 2 . The 
hydroxyl (OH) radical is critical in oxidizing S0 2 to sulfuric 
acid (H 2 S0 4 ) through reactions R } to /? 5 . Because R 4 and R 5 
are relatively fast, sulfuric acid production is limited mainly by 
reactions J x and fl 3 . In this study, the baseline profiles of OH, 
O, and H 2 0 are assumed to be fixed with concentrations taken 
from the photochemical model of Brasseur and Solomon 
[1986]. All the chemical reaction rates correspond to standard 
compilations [DeMore et al. , 1990]. 

Three types of aerosols are considered in this model: 
involatile (or solid) particles which originate from the tropo- 
sphere (we also refer these “tropospheric” CN as Aitken 
particles), pure H 2 S0 4 /H 2 0 particles, and mixed aerosols 
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consisting of Aitken particles in H2SO4/H2O solutions. 
Aitken particles are transported from the troposphere into 
the lower stratosphere and provide the core for mixed 
aerosols either through heterogeneous nucleation or coagu- 
lation between Aitken particles and pure H 2 S0 4 -H 2 0 parti- 
cles. Mixed aerosols interact with H 2 S0 4 vapor by conden- 
sation and evaporation processes. If a mixed aerosol is 
completely evaporated, then its core is released and returned 
to the “Aitken” category. Pure H 2 S0 4 /H 2 0 aerosols can be 
produced only by homogeneous nucleation, when the super- 
saturation of H 2 S0 4 /H 2 0 vapor exceeds a critical value. An 
efficient algorithm to calculate homogeneous and heteroge- 
neous nucleation rates [Zhao, 1993] is used here. Details 
concerning the other microphysical processes and corre- 
sponding numerical algorithms can be found in the works of 
Turco et aL [1979], Toon et al. [1988], and Zhao [1993]. In 
the calculations that follow, each type of particle (Aitken 
nuclei, pure, and mixed particles) is subdivided by size. The 
smallest size treated in the model is about 0.001 fim. The 
total number of particles in our simulations corresponds to 
the sum over the whole size bins. The observed CN corre- 
spond to the total number of particles with sizes greater than 
about 0.01 yarn. To compare with observations, the concen- 
tration of CN in our simulations is by summing all the 
appropriate size bins in the model. 

Simulation Parameters 

Benchmark Simulation 

There are two purposes for investigating the background 
aerosols before simulating transient CN layers. First, we 
need to verify the capability of the model to generate 
realistic vertical concentration profiles of the background 
aerosols and key chemical species. Second, the polar strato- 
spheric model is initialized using the background aerosol 
profiles, which approximate the state of the aerosol layer in 
summer. To achieve a well-characterized background state, 
we first integrate the model to a steady state using 10-year 



S0 2 mixing ratio (ppbv) 

Figure 1. Calculated and observed vertical distributions of 
background S0 2 concentrations. S0 2 concentrations in the 
troposphere are highly variable and cannot be accurately 
predicted with a one-dimensional model. 



Figure 2. Calculated vertical distribution of background 
total aerosol concentration compared to observations [Hof- 
mann, 1988]. 


mean temperatures between 70° and 80°S and a monthly 
mean ozone profile. The tropopause was set at 9 km and was 
kept constant through all of the numerical experiments 
described here. The baseline vertical profiles of S0 2 (Figure 
1) and aerosol concentration (Figure 2) are reasonable com- 
paring with observations. The computed size distributions 
(Figure 3) are similar to those calculated by Turco et al. 
[1979]. The simulated CN concentration is about 5 cm 3 at 
20 km and less than 1 cm -3 at 30 km; these values are close 
to the annual average aerosol concentrations measured near 
the south pole at the same heights. 

Annual Cycle of Temperature 

The observed high correlation between CN concentrations 
and low stratospheric temperatures suggests the importance 
of temperature variation in forcing the anomalous CN. In our 
“polar” simulations the annual cycle of temperature is 
deduced from observed zonal mean temperatures between 
70° and 80°S for 1984 (Figure 4). These data present a typical 
picture of temperature variations in the southern polar 
region; similar temperatures hold in other years. In general, 
temperatures decrease from fall through winter and recover 
in the spring. The lowest temperature in Figure 4 is less than 
190°K, near 18 km in later winter and early spring. An 
approximation to the period of twilight in the polar region is 
superimposed on this temperature field. Near latitude of 
75°S, sunset occurs in early May and sunrise in late August. 
Figure 4 shows clearly that in spring after sunrise, strato- 
spheric temperatures can remain low for a period of several 
weeks. This cycle of temperature and sunlight implies that 
the photochemical production of H 2 S0 4 vapor recovers 
after polar sunrise, while the equilibrium vapor pressure of 
H 2 S0 4 under low-stratospheric temperature is still very 
small. Therefore the acid vapor concentrations can easily 
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Radius (Jim) 

Figure 3. Calculated particle size distributions of the back- 
ground aerosols at altitudes of 17, 21, 25, and 29 km. These 
size distributions compare favorably with measurements 
summarized by Turco et ai [1982]. 

exceed the critical value required for binary homogeneous 
nucleation. Then, new droplets of H 2 S0 4 /H 2 0 solutions can 
be formed and grow to CN sizes by condensation and 
coagulation processes. 

Photolysis of Carbonyl Sulfide (OCS) 

Under ambient stratospheric conditions, sulfate aerosols 
originate from sulfur-bearing gases, especially OCS. Al- 
though one dimensional models are not the best tools for 
simulating OCS transport from the troposphere, the vertical 
stratification can be represented in this way. Indeed, the 
photodecomposition of OCS has been suggested as one 
cause of the CN layer [Rosen and Hofmann , 1983]. OCS 
photolysis is controlled by solar ultraviolet radiation, which 
is particularly sensitive to ozone abundances. Therefore to 
calculate the photodissociation rate of OCS in the southern 
polar region more accurately, monthly average ozone pro- 
files derived from balloon soundings at Syowa (69°S) in 1990 



DAY 

Figure 4. Isotherms of observed zonal belt- mean temper- 
atures between 70° and 80°S latitudes for 1984. The data are 
derived from the National Meteorological Center analysis. 
The shaded area indicates the period of darkness during 
polar winter at a latitude of 75°S. 


were adopted. Stratospheric ozone concentrations decline in 
September, reach minimum values in October, and gradually 
recover in November (Figure 5). These spring ozone “hole” 
profiles have reduced ozone abundances that favor the 
photodecomposition of OCS. 

Baseline CN Simulation 

Using the annual cycles of zonal mean temperature, solar 
illumination, and ozone concentrations derived from mea- 
surements in the southern polar region and the benchmark 
stratosphere as an initial state, a “baseline” simulation of 
CN formation was performed to investigate CN layer forma- 
tion. This baseline model represents the best first-guess 
parameterization for stratospheric CN formation, based on 
existing observations and theories. In the baseline run, 
integration begins at day 80, when the initial state is very 
close to the summer climatology and runs to day 365, when 
the solution essentially returns to the background state. 

The simulated H 2 S0 4 mixing ratio drops dramatically as 
polar sunset advances at around day 120 and remains near 
zero until sunrise the following spring (Figure 6). The loss of 
H 2 S0 4 vapor occurs because the normal balance between 
H 2 S0 4 vapor production and condensation on preexisting 
aerosols cannot be maintained when solar illumination van- 
ishes, and the OH radicals that generate H 2 S0 4 disappear. 
In addition, decreasing temperatures favor H 2 S0 4 loss onto 
the preexisting droplets. Below 35 km, H 2 S0 4 vapor reap- 
pears rapidly after sunrise on day 230; and photochemical 
reactions driven by sunlight reestablish the production of 
H 2 S0 4 . Subsequently, the H 2 S0 4 mixing ratio fluctuates in 
accord with temperature changes. Above about 40 km the 
H 2 S0 4 reduction in midwinter does not occur because of 
higher temperatures there that prevent condensation. Figure 
7 provides several snapshots of the evolution of the vertical 
distribution of H 2 S0 4 concentrations. During the winter 
period the evaporation level for sulfuric acid (i.e., where 
sulfuric acid vapor becomes undersaturated and evaporation 
can occur) rises from 37 to 45 km owing to the wintertime 
trend of temperature decrease in this region. 



Ozone concentration (10 12 cm' 3 ) 

Figure 5. Monthly mean profiles of observed ozone con- 
centrations with unit of 10 * 2 cm' 3 at Syowa Station (69°S, 
40°E) in 1990 [Nankyoku and Kansokittai , 1990]. 
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Figure 8. Variation of total aerosol concentration (number 
per cubic centimeter) with time and altitude in the ‘base- 
line’’ CN simulation. The concentration contours are shown 
with the logarithm with a power of 10. The contour interval 
is 0.5 log units. 


Along with the rapid recovery of H 2 S0 4 just after polar 
sunrise, temperatures in the lower stratosphere are near 
their minimum values (Figure 4). This creates an idea 
situation for particle formation by nucleation, which de- 
pends on H 2 S0 4 concentration and temperature. Indeed, the 
simulated total aerosol concentrations (all particles larger 
than 0.001 /xm) between the tropopause and 35 km increase 
dramatically due to the new particle formation by homoge- 
neous nucleation shortly after sunrise (Figure 8). The aerosol 
concentrations produced by homogeneous nucleation can 
initially exceed 10 3 cm' 3 near 25 km. These newly formed 
sulfate embryos, however, cannot grow to measurable 
“CN” sizes. Only a few particles near 30 km are found 
reaching measurable sizes (greater than 0.01 /xm) and the 
springtime perturbation of CN concentration is hardly no- 
ticeable (Figure 9). The simulated baseline CN profiles are 
not comparable to the observed CN layer (Figure 10). 
Moreover, the aerosol size distributions in the baseline 


experiment show that although a large number of fine 
H 2 S0 4 /H 2 O particles are initially produced by homoge- 
neous nucleation right after sunrise, they do not grow to CN 
sizes because of a shortage of H 2 S0 4 vapor and quick 
coagulation with ambient particles (Figure 11). 

The baseline experiment fails to produce a substantial CN 
layer. The photodissociation rate of OCS is too low to 
produce much H 2 S0 4 vapor in early spring (<10 s ) 
owing to the weak solar illumination. Hence the conversion 
of ambient OCS into H 2 S0 4 cannot be the dominant mech- 
anism in the formation of the CN layer in polar spnng. This 
is true even if the OCS concentration is allowed to accumu- 
late through vertical diffusion over the winter. In fact, the 
polar vortex is highly stable against upward vertical mixing, 
and the vertical profiles of most long-lived tracers suggest 
sustained subsidence in the winter polar vortex, which 
should reduce OCS concentrations at 25-30 km. It follows 
from the baseline simulation that certain mechanisms not 
treated in the baseline experiment must play an important 
role in providing extra sulfur for the formation of new CN. 



H2SO4 concentration (ppbv) 

Figure 7. Vertical profiles of H 2 S0 4 mixing ratio (ppbv) at 
different times is indicated in the figure. The concentrations 
correspond to the baseline CN simulation. Sunrise is at day 
230. 



Figure 9. Variation of CN concentration, i.e., number of 
particles per cubic centimeter with radii exceeding 0.0 1 /xm. 
The CN concentrations are shown as isopleths of constant 
concentration to logarithm. The contour interval is 0.5 log 
units. 
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Aerosol concentration (r > 0.01 p.m) 


of the absence of OH radicals, hence the descending air can 
provide a critical sulfur reservoir for CN formation in the 
spring. The weak CN layer in the baseline experiment is thus 
compensated by a mesospheric sulfur source. 

There are no quantitative observations of the amplitude, 
structure, or duration of the descending motion in the polar 
vortex, but three-dimensional model simulations [Fisher et 
al 1993] provide useful guidelines. The model results show 
that the mean poleward velocity in the mesosphere is about 
40° per month, whereas the subsiding velocity in the upper 
stratosphere of the polar vortex is about 12 km per month in 
July. The descending motion decreases to about 3 km per 
month in the middle stratosphere and the poleward drift is 
also much slower than in the mesosphere. The vortex 
circulation can maintain this structure for up to 6 months. 

On the basis of these model results we postulate a subsid- 
ing velocity field as 



Figure 10. (a) and (b) Vertical distributions of CN concen- 

trations at different times in the baseline model simulation, 
compared to observations [Hofmann and Rosen , 1985]. 


Role of Subsiding Motion in the Polar Vortex 

A Lagrangian visualization using a three-dimensional 
model simulation provides a detailed picture of the airflow in 
the middle atmosphere [Fisher et aL , 1993]. Air from the 
mesosphere executes an inward and downward spiral into 
the polar vortex. The numerical simulation offers an expla- 
nation of HALOE observations of low mixing ratios of 
methane in the polar stratosphere [Russell et al 1993]. The 
subsiding motion in the polar vortex transports mesospheric 
air with low mixing ratios of methane down into the strato- 
sphere. Similarly, the subsiding motion would transport a 
higher mixing ratio of sulfur dioxide into the polar strato- 
sphere from the mesosphere. A high mixing ratio of S0 2 in 
the mesosphere is likely from the photolysis of H 2 S0 4 at 
those altitudes [Turco et al . , 1979]. The S0 2 transported 
downward would not be destroyed in the polar night because 



t - 220 
180 


120 < / <310. 


( 1 ) 


Here t is time (day), p T is the air density at the top layer of 
the model, vr 0 = 0.5 cm/s, and a = 0.4. This idealized 
formulation of vertical descent qualitatively agrees with the 
results of the three-dimensional simulation. 

The effect of the corresponding poleward airflow is treated 
consistently by including an additional term in the one- 
dimensional model to ensure mass continuity. First, note 
that the horizontal advection is related to the vertical veloc- 
ity divergence though the mass continuity equation. More- 
over, the continuity equation can be written for any tracer in 
the form 


dC d(\rC) , C\ 

— + • (v H C) + — = V • pKV — ^ P - L, (2) 

dt dz \p J 



1.0E-3 1.0E-2 1.0E-I 1.0E+Q 1.0E+I 

Radius (|im) 


Figure 11. Aerosol size distributions at different times in 
the baseline model simulation. 
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If zonal symmetry is assumed, the horizontal advection term 
can be approximated by differencing from the pole over the 
radius / of the vortex as 


V/r (v H C) 


V hC a 

T~ 


(3) 


Here the subscript a indicates values in ambient air outside 
the vortex at the same altitude, and v h is the average 
horizontal poleward velocity. Similarly, the mass continuity 
equation can be written as 


v H 1 d(wp) 

+ = 0 . 

/ p dz 


(4) 


Then, in a one-dimensional tracer representation we have 


dC d(wC) 

h 

dt dz 


= C a 


1 B(wp) d 


d C 


dZ 


dz 


dZ \p 


+ -P*77 7 


(5) 


The first term on the right-hand side of (5) can be viewed as 
the effect of horizontal mass convergence associated with 
net vertical subsidence. Whereas the vertical velocity used 
here increases with altitude, the downward mass flux actu- 
ally decreases with altitude, and the horizontal convergence 
effectively brings air into the vortex from outside. Both the 
vertical motion and the horizontal convergence decrease in 
the lower stratosphere, and there is little accumulation of 
mass in the stratospheric column. 

As a variation of a baseline experiment, a S0 2 mixing ratio 
of 0. 1 parts per billion by volume (ppbv) is assumed at the 
upper boundary (60 km). Similarly, the mixing ratio of 
H 2 S0 4 in air descending from mesosphere is set to 0.1 ppbv, 
and the concentrations of preexisting aerosols are set to zero. 
These concentrations for S0 2 and H 2 S0 4 are taken from 
photochemical simulations of the stratospheric sulfur cycle 
[e.g., Turco et al . , 1979, 1982]; however, the actual concentra- 
tions are uncertain. As will be shown, the results are not 
sensitive to the H 2 S0 4 concentration but are sensitive to S0 2 . 



Aerosol concentration (r > 0.01 jim) 


Figure 12. Vertical distributions of CN concentrations at 
different times in the experiment including a subsiding mo- 
tion from the mesosphere. Compare with data shown in 
Figure 10. 



Aerosol concentration (r > 0. 15 gm) 


Figure 13. The vertical distribution of large aerosol parti- 
cles with r > 0.15 p,m, versus time in the simulation with 
subsidence. 


After including such a descending motion and associated 
fluxes at the upper boundary, the model generates dense CN 
layers during spring after sunrise that agree with observa- 
tional data (Figure 12). The layer appears at higher altitudes 
in early spring and shifts downward as time progresses. The 
maximum concentration of CN particles is situated near 29 
km in September and 24 km in October. The background 
particles in this case also have lower concentrations during 
winter and spring because of the downward flow, which 
brings particle-free air into the upper and middle strato- 
sphere. The number of large particles is reduced significantly 
(Figure 13). Hence subsidence is also consistent with the 
aerosol hole observed by balloon-borne and satellite sound- 
ings [Hofmann and Rosen , 1983; McCormick et al., 1983]. 

The vertical flux of S0 2 at upper boundary' is crucial to 
forming the CN layer. During fall the photolysis of OCS 
decreases and the S0 2 mixing ratio above 20 km decreases 
through reaction with OH. The photochemical source of S0 2 
remains insignificant throughout the winter. Moreover, the 
production of S0 2 from OCS is too slow to generate a CN 
layer in the springtime. Once the downward transport of S0 2 
from the mesosphere into the stratospheric vortex is in- 
cluded, however, the S0 2 mixing ratio can be maintained at 
an elevated level, as shown in Figure 14. 

Variations in the H 2 S0 4 mixing ratios for the case with 
subsidence are similar to those in Figure 6, except the mixing 
ratio of H 2 S0 4 has a higher value near 30 km in winter and 
early spring due to the subsidence motion; however, the 
total aerosol and CN concentrations (Figures 15 and 16, 
respectively) are considerably different from those in Fig- 
ures 8 and 9. The perturbations in the aerosol concentrations 
following sunrise are far more intense and persistent. Note 
particularly the strong enhancement in CN concentrations 
from 18 to 30 km during springtime in Figure 16. By 
comparing Figures 15 and 16, it becomes evident that the 
enhancement in CN concentrations occurs mainly above 18 
km, although many small droplets are generated below this 
level (these very small particles are not counted as CN). At 
the lower altitudes the ambient aerosols capture all of the 
small droplets by coagulation before they grow to CN sizes 
by condensation. The vertical distribution of the CN con- 
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SO2 concentration (ppbv) 


Figure 14. Vertical profiles of S0 2 mixing ratio (ppbv) at 
different times of the year for the case with subsidence. 

centration near the evaporation level at 35 km (Figure 15) is 
quite sensitive to local temperature variation. A clear corre- 
lation between fluctuations in the temperature and simulated 
CN is seen during springtime. 

The simulated aerosol size distributions for the case with 
subsidence reveal the detailed evolution of these aerosols 
(Figure 17). The aerosol size distribution in winter is similar 
to the ambient size distribution, as might be expected, 
although the size spectrum shifts to smaller-sized concentra- 
tions as some of the bigger particles fall out. Ten days after 
sunrise, however, a large number of very small H 2 S0 4 /H 2 0 
particles are generated by homogeneous nucieation, and the 
size distribution is bimodal. Unlike the baseline case, the 
aerosol size spectra for days 50 and 80 (after sunrise) 
indicate that the small particles have grown to larger sizes, 
and the typical particle radius on day 80 is at 0.01 /xm. 

Figure 18 shows the complete evolution of the CN distri- 
bution from sunset through the end of the year at altitudes 
from 10 to 60 km. Clearly, a significant enhancement in the 
CN concentration occurs during spring (September, Octo- 
ber, and November) after sunrise. The peak in the CN layer 



Time (day) 


Figure 15. Variation^ total aerosol concentration (num- 
ber per cubic centimeter) with time and altitude in the 
experiment with a subsiding motion. The concentration 
contours are shown with the logarithm of the concentrations. 
The contour interval is 0.5 log units. 



Figure 16. Variation of CN concentration (number per 
cubic centimeter) with time and altitude in the experiment 
with a subsiding motion. The CN concentrations are shown 
as isopleths of constant concentration to logarithm. The 
contour interval is 0.5 log units. 


moves from the highest altitude early in spring down to 23 km 
in October. The maximum concentration of CN particles 
reaches 40 cm' 3 , with the main layer located between 18 and 
37 km. Below 15 km, CN-sized particles are lost by coagulating 
with the larger number of ambient aerosols. Above 37 km, 
higher temperatures cause the volatile CN droplets to evapo- 
rate. The simulated CN layer is consistent with measurements 
of the intensity, altitude range, and time of occurrence. 

In the simulations, subsiding air carries higher mixing 
ratios of S0 2 from the mesosphere into the stratospheric 
polar vortex thereby creating a CN layer, which seems to 
provide a reasonable explanation for CN layer formation and 
the springtime aerosol hole. The case with subsidence also 
reproduces the observed major features and evolution of the 
polar CN layer. We conducted a series of sensitivity exper- 
iments with different subsidence velocities and mesospheric 
S0 2 mixing ratios. The intensity of the CN layer increases 
with higher S0 2 mixing ratios and rates of descending air. It 



Radius (|xm) 

Figure 17. Aerosol size distributions at different times in 
the experiment with a subsiding motion. 
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Figure 18. Vertical and temporal distribution of CN con- 
centration from days 120 to 360 at altitudes between 15 and 
59 in the simulation with subsidence. 


is found that with an S0 2 mixing ratio of 0.1 ppbv or larger, 
a parameter a in (1) of 0.4 ~ 0.5 and a velocity w 0 in the 
range 0.4 ~ 0.5 cm s' 1 , the generated CN layers are similar 
to those observed. According to the observational data from 
Toon et al. [1989] the vertical velocity near 18 km is about 
-85 m d' 1 (-0.1 cm s' 1 ). The range of maximum subsi- 
dence rates near 18 km used in the model is between 0.02 and 
0.1 cm s" 1 which is close to the observations. 

We also run the case with no S0 2 but 0.1 ppbv H 2 S0 4 
mixing ratio at the model upper boundary. We found that the 
simulated CN layer is very weak; the peak concentration is 
less than 10 cm' 3 in this simulation. Moreover, the strength 
of the layer is not sensitive to increasing the H 2 S0 4 concen- 
tration at the upper boundary. The lack of sensitivity occurs 
because H 2 S0 4 vapor can condense on preexisting particles 
as it is transported from high altitudes to the polar strato- 
sphere. The short lifetime of H 2 S0 4 vapor prevents its 
accumulation in the region where CN events occur. In 
contrast, the S0 2 is not destroyed in the polar night and 
therefore can serve as a sulfur reservoir for the CN forma- 
tion in springtime. 


Modified Photochemistry 

To test the importance of oxidation of OCS by OH radicals 
originally proposed by Oppenheimer [1987] and other related 
chemical processes, the following photochemical reactions 
involving OCS were considered. 

(J { ) OCS + hv-+ S + CO 

( R 6 ) O + OCS -* SO + CO, K"[0] ~ 10 - U (cm 3 s' 1 ) 
(/? 7 ) OH + OCS -*SH + CO 2 , K'iohi ~ !0 -15 (cm 3 s' 1 ) 
(R h ) CIO + OCS -» products, ^[cio] < 10 16 (cm 3 s ') 
(/? 9 ) N0 3 + OCS -* products, AT NOj ] < 10 15 (cm 3 s ’) 
(/?,„) Cl + OCS — SCI + CO, /sT fcI ] < 10~ 16 (cm 3 s' 1 ) 

( 6 ) 


The first two reactions are included in the baseline model. As 
we noted earlier, the photolysis rate of OCS (/ t ) is too small 
in polar spring to provide an adequate sulfur source for the 
observed CN layer. An OCS photodecomposition rate, /ocs, 
of roughly 10"' s' 1 is necessary to create a realistic CN 
layer. At 25 km, the reaction rate with atom oxygen, / 6 = 
/? 6 [0] is about 10" 8 (if [O] ~ 10 6 cm -3 [Brasseur and 
Solomon, 1986]. This is obviously too slow. The reaction of 
OH with OCS (Rj) suggested by Oppenheimer [1987], with 
an upper limit concentration for OH of — 10 7 cm' 3 , is again 
too small (/ 0 h ~ 10 ' 8 s' 1 ). Measurements show that CIO 
mixing ratios in the lower polar stratosphere during spring 
can exceed 1 ppbv, more than 100 times greater than normal 
background stratospheric concentrations [de Zafra et al., 
1987; Solomon et al ., 1987]. Thus the OCS lifetime against 
decomposition by reaction R% could become as large as 
about 10 “ 7 s' 1 , which is an order of magnitude greater than 
reaction / 6 . This chemical reaction may contribute to the 
conversion of OCS to H 2 S0 4 vapor. /?9 and R\$ would be 
more than an order of magnitude less important than /? 8 , 
because the concentrations of N0 3 and Cl are so low. 

A CIO concentration field is constructed according to 
observational data from de Zafra et al. [1987] and Solomon 
et al. [1987] as 

j , 230 < t < 270 (7) 

where 


[CIO] = c(z) cos 


c(z) = 


exp 

(z-20) 4 ' 

2 

exp 

i. 



z < 20 km 


z > 20 km 


( 8 ) 


Here t is time (day) and [CIO] is CIO mixing ratio (ppbv). 

When we incorporate reaction R% in the baseline model 
with the CIO concentrations shown in Figure 19 and the 
upper limit (room temperature) rate coefficient of 10 16 cm 
s' , CN layers are formed between 15 and 30 km with peak 
concentrations of more than 100 cm 3 (Figure 20). This 
surplus of CN reflects particle formation, condensation, and 
coagulation in the presence of abundant H 2 S0 4 vapor gen- 
erated by the decomposition of OCS. 

Among the photochemical processes listed above, reac- 
tion is the most likely to provide sufficient H 2 S0 4 vapor 
for CN formation in polar spring. The reaction rate coeffi- 
cient, however, is probably much smaller than the room 
temperature upper limit [DeMore et al., 1990] used in the 
simulation. At much lower stratospheric temperatures than 
the room temperature the rate coefficient may be 1 to 2 
orders of magnitude smaller. Further experimental studies of 
this mechanism are needed to eliminate this possibility. 


Effect of Micrometeorites 

To study the impact on CN layer formation of an addi- 
tional particle source associated with micrometeorites, we 
also did some experiments which include a constant flux of 
small particles at the upper boundary according to the results 
from Hunten et al. [1980] and Turco et al. [1980]. These fine 
meteoritic aerosols slowly migrate downward by diffusion in 
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Time (day) 


Vertical profile at day 240 



CIO mixing ratio (ppbv) 

Figure 19. (a)-(c) CIO mixing ratios used in the numerical experiment with reaction R % set to its upper 
limit value. The enhanced CIO concentrations were constructed using equation (7) referring to the 
measurements from de Zafra et al. [1992]. 


the one-dimensional simulation. Figure 21 shows the steady 
state aerosol number size distributions obtained by integrat- 
ing the model 300 days under conditions corresponding to 
the background state with a meteoritic aerosol source at the 
upper boundary. Even though there are a larger number of 
meteoritic dust particles throughout the upper stratosphere, 
these do not influence CN concentrations. The aerosol size 
spectra (Figure 21) show a micrometeorite mode below a 
0.01-^tm radius, with a typical radius of 0.003 ^m. The initial 
size of the meteoric dust particles, although quite uncertain, 
are very unlikely to be larger than 0.01 fxm [Hunten et al ., 
1980]. Accordingly, micrometeoritic particles themselves 
are unlikely to change the CN (r > 0.01 ^im) concentration, 
or to form a distinct CN layer. 

The micrometeoritic particles were incorporated in the 
CN simulation as an initial distribution (Figure 21) in equi- 


librium with a fixed flux at the upper boundary. Other 
conditions were similar to those in the baseline experiment. 
The resultant CN concentrations, including the vertical 
distributions and temporal evolution, are almost the same as 
in the baseline model. When micrometeorites are coupled 
with either subsiding motion or reaction R s in the model, 
minor variations in CN concentrations occurred when het- 
erogeneous nucleation was assumed to have a contact angle 
of 50° (moderate nuclei). The overall features of the CN 
layers remained unchanged. Clearly, the influence of mi- 
crometeorites is a secondary factor in CN layer formation. 

We also considered an extreme case by assuming (1) all 
the micrometeorites serve as very efficient condensation 
nuclei for sulfate aerosols and are nucleated in a very short 
time once H 2 S0 4 vapor is supersaturated and (2) homoge- 
neous nucleation is switched off so that no new particles are 
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generated. When we incorporated conditions 1 and 2 into the 
model with otherwise the same conditions as in baseline 
experiment, the model results are almost the same as those 
in baseline experiment. A CN layer still does not appear. 
Even if the micrometeorites are assumed to be excellent 
nuclei for sulfuric acid and water solutions, they did not 
enhance CN formation in this experiment. 

When conditions l and 2 are combined with vortex subsi- 
dence and sulfur flux at the upper boundary of the model, a 
CN layer is produced similar to that shown in Figure 12. 
However, the scenario of CN formation is distinctly different 
from that in the similar case without meteorites. In the case 
illustrated in Figure 12 the new sulfuric acid/water droplets 
are produced by homogeneous nucleation. In contrast, in 
this case, it is the heterogeneous nucleation of sulfuric 
acid/water solutions on these preexisting micrometeorites 
that is the key mechanism to initialize the CN formation. As 
we can see in Figure 21, the concentrations of micromete- 
orites can be more than 100 cm 3 above 30 km. These 
micrometeorites are heterogeneously nucleated and grow to 
CN sizes (>0.01 n m) by condensation as they are trans- 
ported to the lower stratosphere by subsidence. 

According to the results from numerous numerical exper- 
iments we found that if sufficient concentrations of wettable 
micrometeorites exist, we can not rule out the possibility of 
CN layer formation through heterogeneous nucleation of 
sulfuric acid and water solutions on these meteorites. How- 
ever, the presence of sufficient H 2 S0 4 vapor concentrations 
is the most crucial factor controlling the occurrence of CN 
events. With the subsidence supplying the extra sulfur from 
the mesosphere, a CN layer can be generated either through 
homogenous or heterogeneous nucleation, or by both. 

Influence of Polar Stratospheric Clouds (PSCs) 

It has been proposed that the largest sulfate particles in the 
stratosphere may freeze and serve as condensation nuclei for 
PSCs. The removal of aerosols from the polar vortex may be 
related to PSC formation in winter and early spring when 
temperatures drop below about 193°-195°K (refer to the 



Aerosol concentration (r > 0.01 |im) 

Figure 20. Vertical distributions of CN concentrations at 
different times in the numerical experiment with reaction /? 8 
set to its upper limit value and elevated CIO concentrations 
(see Figure 19). 



Radius (\im) 

Figure 21. Simulated aerosol size distributions for back- 
ground conditions including meteoric smoke particles [after 
Hunten et al . , 1980]. 


temperature variation in Figure 4). The larger PSC particles 
fall to lower altitudes more quickly than sulfate aerosols and 
thereby serve as a sink. With the largest sulfate aerosols 
removed, H 2 S0 4 vapor loss by condensation on particle 
surfaces is reduced, which may favor new particle formation 
by nucleation. 

To test the influence of PSCs on CN formation, it was 
assumed that whenever the temperature falls below 195°K, 
large sulfate particles (r > 2.5 /u.m) are scavenged by PSCs. 
The results suggest that the loss of aerosol surface area due 
to PSC sedimentation during winter is not a critical mecha- 
nism for CN formation. That is, the baseline experiment 
does not generate a CN layer even if this process is included; 
however, when it is combined with vortex subsidence or the 
CIO reaction, the predicted CN layers are slightly intensi- 
fied. 


Summary and Discussion 

The major conclusions derived from numerous model 
experiments regarding polar CN layer formation are summa- 
rized in Table 2. We found that the variations in photochem- 
ical reactions corresponding to the conditions of polar twi- 
light and low temperature control the first appearance of new 
particle formation. Specifically, OH radicals play an impor- 
tant role in producing sulfuric acid vapor from S0 2 , and low 
temperatures are crucial for homogeneous nucleation, which 
generates large numbers of small H 2 S0 4 /H 2 0 particles. 
However, If OCS photolysis is the only mechanism consid- 
ered, the newly formed particles cannot grow to observable 
CN sizes because of the shortage of H 2 S0 4 vapor. Not only 
is the photodissociation rate of OCS too small under the 
weak solar illumination at large zenith angles in early spring 
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Table 2. Summary of the Importance of Various Mechanisms for CN Production 


Mechanism 

Contribution to CN Layer Formation 

Low temperatures 
OCS photolysis 
Vortex dynamics 
Oxidation of OCS by OH 
Reaction of OCS with CIO 
Micrometeorite influx 
PSC scavenging 

essential (initiates homogeneous nucleation) 

secondary (supplies some H 2 S0 4 vapor) 

essential (provides high S0 2 mixing ratios through subsidence) 

unimportant (unless OH concentrations exceed 10 7 cm" 3 ) 

possibly unimportant (unless the reaction has upper limit value and a small action energy) 
possibly important (if sufficient concentrations exist and are active as condensation nuclei) 
minor (except perhaps in the Antarctic vortex where substantial PSC-induced sulfate sedimentation may occur) 

CN, condensation nuclei; OCS, carbonyl sulfide; PSC, polar stratospheric cloud. 


but also the OCS concentrations would be suppressed by 
subsidence in the polar vortex over the winter season. 

Air descending from the mesosphere in the polar vortex 
carries higher S0 2 mixing ratios into the stratosphere and 
provides an extra sulfur source for the formation of a CN 
layer. In a one-dimensional model, realistic CN layers are 
produced when the vortex dynamical process is included. 
The major features of the CN layer observed by balloon- 
borne instruments are accurately depicted in the “subsi- 
dence” experiment. For instance, the layer appears at higher 
altitudes in early austral spring and shifts downward to near 
25 km by October. The layers are confined to altitudes 
between 18 and 37 km. Below 18 km, the ambient particles 
absorb sulfuric acid vapor fast enough to prevent new CN 
production. Above 37 km the local temperatures are high 
enough to prevent H 2 S0 4 vapor from nucleating. The sub- 
sidence experiment also explains the small concentrations of 
optically active aerosols during winter and spring resulting 
from the intrusion of aerosol-free air from high altitudes. 
Polar subsidence and photochemistry of S0 2 seems to be the 
most plausible explanation for CN layer formation in polar 
spring. Our calculations are based on the parameters such as 
temperature, subsidence, and solar illumination for the Ant- 
arctic stratosphere, but these conditions are similar over 
Arctic. We believe CN events in both hemispheres have 
similar formation mechanisms. 

By analyzing all of the reactions involving stratospheric 
OCS that might initiate sulfate aerosol production in spring, 
we have found that none of these reactions are likely to 
provide sufficient H 2 S0 4 vapor for CN layer formation and 
persistence. The photochemical production of S0 2 by the 
reaction between OCS and CIO seems to be a secondary 
process at best. The OCS + CIO reaction, although quite 
uncertain in the absence of sufficient laboratory data, would, 
at its upper limit room temperature value, be barely ade- 
quate to generate a weak CN layer. 

Other sensitivity tests were also conducted. We conclude, 
for example, that micrometeorites may cause minor varia- 
tions in CN concentrations. Heterogeneous nucleation on 
micrometeorites could be important for CN formation if the 
micrometeorites are good condensation nuclei for sulfuric 
acid and water solutions and have sufficient concentrations 
in the stratosphere. Homogeneous nucleation in most cases 
appears to be necessary to explain the increase in CN 
concentrations. Indeed, the total sulfate mass in the CN 
layer implies such a large H 2 S0 4 vapor production rate that 
homogeneous nucleation at polar temperature seems un- 
avoidable. Reductions in the background sulfate aerosol 
surface area caused by PSC condensation and sedimentation 
can help CN layer formation according to our model simu- 


lations. However, PCS effects cause minor variations in CN 
layer properties and cannot explain the appearance of the 
layer or its overall characteristics. 
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